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Abstract. Relativistic supernovae have been discovered until recently only through their association with long duration 
Gamma Ray Bursts (GRB). As the ejecta mass is negligible in comparison to the swept up mass, the blastwaves of such 
explosions are well described by the Blandford-McKee (in the ultra relativistic regime) and Sedov-Taylor (in the non- 
relativistic regime) solutions during their afterglows. However, the recent discovery of the relativistic supernova SN 2009bb, 
without a detected GRB, has indicated the possibility of highly baryon loaded mildly relativistic outflows which remains in 
nearly free expansion phase during the radio afterglow. In this work, we consider the dynamics and emission from a massive, 
relativistic shell, launched by a Central Engine Driven EXplosion (CEDEX), decelerating adiabatically due to its collision 
with the pre-explosion circumstellar wind profile of the progenitor. We show that this model explains the observed radio 
evolution of the prototypical SN 2009bb and demonstrate that SN 2009bb had a highly baryon loaded, mildly relativistic 
outflow. 
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INTRODUCTION 

Gamma-ray Bursts (GRBs) have long been recognized to require ultra relativistic bulk motion of matter particles 
(see Piran [1,2] for reviews). GRB afterglows are generated from the emission by relativistic shocks that result from 
slowing down of a relativistic shell by the the medium surrounding the progenitor star that exploded. In core collapse 
supernovae similar interaction of stellar material (ejecta) from an exploding star with the circumstellar matter (CSM) 
results in non-relativistic shocks. 

Fluid dynamics of ultra-relativistic spherical blast waves mediated by strong shocks has been treated by Blandford 
and McKee [3, 4]. They found a similarity solution of an explosion of a fixed amount of energy in a uniform 
medium. On the other hand. Chevalier [5], Nadezhin [6] described the initial nearly free expansion of a non-relativistic 
supernova blast wave, interacting with the surrounding circumstellar medium. Once the blast wave sweeps up more 
CSM material than its own rest mass, the self-similar solutions of non-relativistic blast waves are described in the 
Newtonian regime by the Sedov [7] von Neumann [8] Taylor [9] solution. 

In this presentation we provide an analytic solution (see ApJ [10] for details) of the standard model of relativistic 
hydrodynamics [see e.g. 1, 11] for an adiabatic blastwave. Here, the exploding shell decelerates due to inelastic 
collision with an external medium. The solution provided here is for an arbitrary Lorentz factor of the expanding 
supernova shell. This solution which can handle a trans-relativistic outflow is motivated by the discovery of SN 
2009bb, a type Ibc supernova without a detected GRB which shows clear evidence of a mildly relativistic outflow 
powered by a central engine [12]. SN 2009bb-like objects (Central Engine Driven Explosions, hereafter CEDEX [10]) 
differ in another significant way from classical GRBs: our work shows that they are highly baryon loaded explosions 
with non-negligible ejecta masses. The new analytic blastwave solution here therefore generalizes the Blandford and 
McKee [3] result, in particular their impulsive, adiabatic blast wave in a wind-like p r^^ CSM. 



RELATIVISTIC BLASTWAVE SOLUTION 

We use the simple collisional model described by Piran [1], Chiang and Dermer [11] where the relativistic ejecta 
forms a shell which decelerates through infinitesimal inelastic collisions with the circumstellar wind profile. The 
initial conditions are characterized by the rest frame mass Mq of the shell launched by a CEDEX and its initial Lorentz 



factor 7o- 

The shell slows down by collision with the circumstellar matter. The swept up circumstellar matter is given by m{R). 
Conservation of energy and momentum give us [see 11, 1], 
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respectively, where dE is the kinetic energy converted into thermal energy, that is the energy in random motions as 
opposed to bulk flow, by the infinitesimal collision. 

For a circumstellar medium set up by a steady wind, where we expect a profile with p r^^, we have m{R) = AR 
where A is the mass swept up by a sphere per unit radial distance. A = M /vwind can be set up by a steady mass loss 
rate of M with a velocity of v„ind from the pre-explosion CEDEX progenitor, possibly a Wolf Rayet star Integrating 
the right hand side and solving for 7 > 1 this equation we have 
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which gives the evolution of 7 as a function of R. The amount of kinetic energy converted into thermal energy when 
the shell reaches a particular R can be obtained by integrating Equation (1) after substituting for 7 from Equation (2) 
and dm = AdR, to get 
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E^c"- y-MQ -AR+ ^JM^ + IAjqRMq +A^R^j . (3) 

The evolution of R and 7 can be compared with observations once we have the time in the observer's frame that 
corresponds to the computed R and 7. For emission along the line of sight from a blastwave with a constant 7 the 
commonly used expression [13] is tots ~ R/i^Y^c). However, Sari [14] has pointed out that for a decelerating ultra- 
relativistic blastwave the correct is given by the differential equation dto^s — dR/{2'f-c). We substitute 7 from 
Equation (2) and integrate both sides to get the exact expression 
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Note that, this reduces to the [13] expression only in the case of nearly free expansion and deviates as the shell 
decelerates. In the rest of the work we use t to indicate the time f„/,i in the observer's frame. Inverting this equation and 
choosing the physically relevant growing branch, gives us the analytical time evolution of the line of sight blastwave 
radius, as 

1 



R = ^ X i^-Mo + lAcjQt + ^ McM^tjl + (Mo - 2Ac7of )2 j , (5) 

in the ultra-relativistic regime. This can now be substituted into Equations (2 and 3) to get the time evolution of 
the Lorentz factor 7 and the thermal energy E [10]. This completes the solution for the blastwave time evolution, 
parametrized by the values for 70, Mo and A. 



BLASTWAVE ENERGETICS 

Chakraborti and Ray [10] use the blastwave solution developed in the previous section to predict the radio evolution of 
a CEDEX with a relativistic blastwave slowing down due to circumstellar interaction. For the prototypical SN 2009bb, 
the blastwave was only mildly relativistic at the time of the observed radio afterglow. In the absence of a significant 
relativistic beaming, the observer would receive emission from the entire shell of apparent lateral extent at a 
time tgf,^ given by dt^f,^ = dRi^t/(j5Yc)- Integrating term by term gives us the time evolution of the lateral radius as 

Riar = cy^7Q — If — O (t^y The thermal energy available when the shell has moved out to a radius R is given exactly 

by Equation (3), however it is again convenient to look at its Taylor expansion, E =Ac^{yo — 1)R — {R'^)- 
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FIGURE 1. A: Evolution of the blast wave radius Riat, determined from SSA fit to observed radio spectrum, as a function of 
the observer frame time fo^,. The evolution is consistent with nearly free expansion and observations require Mq > lO^^'^M,-:,. B: 
Magnetic field as a function of blast wave radius, as determined from SSA fits. Blue dots represent size and magnetic field of radio 
supemovae from Chevalier [15] at peak radio luminosity. Red crosses (with 3(T error-bars) give the size and magnetic field of SN 
2009bb at different epochs, from spectral SSA fits. Red line gives the best B oc R^^ (Equation 7) fit. 



We consider a electron distribution with an energy spectrum NqE^^cIE, which we assume for simplicity to be 
extending from j„,meC^ to infinity, filling a fraction / of the spherical volume of radius R. If a fraction £e = Ee/E of 
the available thermal energy goes into accelerating these electrons, then for the leading order expansion of E mR the 
normalization of the electron distribution is given by 

_ 3Ac^£,(7o-l)(7,„m,c^)^ 
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for /:> = 3, as inferred from the optically thin radio spectrum of SN 2009bb [12]. 

We consider a magnetic field of characteristic strength B filling the same volume. If a fraction Eg = Eg/E goes into 
the magnetic energy density, then the characteristic magnetic field is given by 



This explains the observed B °^ R^ behavior (Figure 1) seen in SN 2009bb. These observations strengthens the case 
for a nearly constant eg. 

Note that the highest energy to which a cosmic ray proton can be accelerated is determined by the B R product 
[16, 17]. We have argued elsewhere [18] that the mildly relativistic CEDEX are ideal for accelerating nuclei to the 
highest energies to explain the post GZK cosmic rays. 



EXTRACTING BLAST WAVE PARAMETERS FROM RADIO OBSERVATIONS 

The inverse problem is that of determining the initial bulk Lorentz factor specified by 70, progenitor mass loss rate 
given by A or M and the initial ejecta mass Mo, from the radio observations. The bulk Lorentz factor may be determined 
from the radio observations using to get the simplified expression for 70 as [10], 
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The result is insensitive to the equipartition parameter a = £e/£B and filling fraction /. This may be used to reliably 
determine the initial bulk Lorentz factor of a radio detected CEDEX in the mildly relativistic, nearly free expansion 
phase (like SN 2009bb). 

A simplified expression for A = M /v^^ind, the circumstellar density profile, set up by the mass loss from the 
progenitor has been derived by Chakraborti and Ray [10]. This gives us the approximate expression for the mass 



loss rate as 
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This approximate expression indicates the dependence of the inferred mass loss rate on the observational parameters, 
and makes an error of only < 10% in case of SN 2009bb, when compared to our exact expression. Note that, this 
expression has similar scaling relations as Equation (23) of Chevalier and Fransson [19]. Hence, the mass loss rate of 
SN 2009bb as determined using that equation by Soderberg et al. [12] remains approximately correct. 

The initial ejecta rest mass Mq cannot be estimated from radio observations in the nearly free expansion phase. It can 
only be determined when the CEDEX ejecta slows down sufficiently due to interaction with the circumstellar matter 
(Figure 1). Thereafter, the initial ejecta mass can be obtained using the timescale of slowdown fj^^ and the already 
determined A and Yq. Nearly free expansion for a particular period of time, can only put lower limits on the ejecta 
mass, as shown in this work. 



DISCUSSIONS 

Chakraborti and Ray [10] provide a solution of the relativistic hydrodynamics equations which is uniquely tuned to the 
CEDEX class of objects hke SN 2009bb. Because of the non-negligible initial ejecta mass of such a CEDEX, objects 
like this would persist in the free expansion phase for quite a long time into their afterglow. Sweeping up a mass equal 
to that of the original ejecta would take considerable time, unless the mass loss scale in its progenitor was very intense, 
i.e. it had a large A. We refer the reader to Chakraborti and Ray [10] for a comparison of the CEDEX solution to 
the intermediate Blandford-McKee like solution and a final Snowplough phase. We also invert the dependence of the 
observed parameters on Yo and M in terms of the peak frequency, peak time and peak fluxes to interpret the parameters 
of the explosion. 
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